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Abstract 

All possible symmetry-determined nonlinear normal modes (also called by simple periodic orbits, one-mode solutions etc.) 
in both hard and soft Fermi-Pasta-Ulam-/3 chains are discussed. A general method for studying their stability in the thermo- 
dynamic limit, as well as its application for each of the above nonlinear normal modes are presented. 
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I. INTRODUCTION 

This paper is devoted to studying the stability of dy- 
namical objects which are called by very different terms 
such as one- mode solutions (OMSs) [Ij, simple peri- 
odic orbits (SPOs) [sl, low-dimensional solutions [1], one- 
dimensional bushes [1, Q etc. Below we refer to them as 
nonlinear normal modes (NNMs). Let us comment on 
this terminology. 

The concept of similar nonlinear normal modes was 
developed by Rosenberg many years ago 12] . Each NNM 
represents a periodic vibrational regime in the conser- 
vative A'^-particle mechanical system for which the dis- 
placement Xi{t) of every particle is proportional to the 
displacement of an arbitrary chosen particle, say, the first 
particle [a:i(t)] at any instant t: 

Xi{t) = Ci ■ Xl{t), (1) 

where {ci = 1, C2, C3, . . . , cn} are constant 
coefficients. [27} 

Note that convenient linear normal modes (LNMs) also 
satisfy Eq. ^ since, for any such mode, one can write 



Xi{t) — Oi ■ sin(a;t + (po), i — 



(2) 



where are constant amplitudes of individual particles, 
while w and (j)o are the frequency and initial phase of the 
considered mode. 

As a rule, NNMs can exist in the mechanical systems 
with rather specific interparticle interactions, for exam- 
ple, in systems whose potential energy represents a ho- 
mogeneous function with respect to all its arguments. 
However, in some cases, the existence of NNMs is caused 
by certain symmetry-related reasons. We refer to such 
dynamical objects as symmetry-determined NNMs. 

In Q, we have found all symmetry-determined NNMs 
in all iV-particle mechanical systems with any of 230 
space groups. This proved to be possible due to the 
group-theoretical methods developed in [9l-fTll| for con- 
structing bushes of vibrational modes. 



At this point, it is worth to comment on the concept 
of bushes of modes introduced in [9]_(the theory of these 
dynamical objects can be found in fill Fl^). 

In rigorous mathematical sense, they represent 
symmetry-determined invariant manifolds decomposed 
into the basis vectors of irreducible representations of the 
symmetry group characterizing the considered mechani- 
cal system ( "parent" group) . 

Because of the specific subject of the present paper, it 
is sufficient to consider only bushes of vibrational modes 
in nonlinear monoatomic chains. Such bushes have been 
discussed in d, [1] . Let us reproduce here some ideas and 
results from these papers. 

Every bush B[G'] describes a certain vibrational regime, 
and some specific pattern of instantaneous displacements 
of all the particles of the mechanical system corresponds 
to it. In turn, this pattern is characterized by a symmetry 
group G (in particular, such group can be trivial) which 
is a subgroup of the symmetry group Gq of the mechanical 
system in its equilibrium state. 

For example, let us consider the two-dimensional bush 
B[a^,i] in the monoatomic chain with periodic bound- 
ary conditions whose displacement pattern X(t) can be 
written as follows 



Xit) ^{x,{t),X2it),-X2{t),~Xi{t) 
-X2{i),-Xi{t) I ...}. 



Xi{t),X2{t), 



(3) 
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This pattern is determined by two time-dependent func- 
tions Xi{t), X2{t), and the corresponding vibrational state 
of the A^-particle chain is fully described by displace- 
ments inside the primitive cell, which is four time larger 
than that of the equilibrium state. 

We will refer to the ratio of the primitive cell size of 
the vibrational state to that of the equilibrium state as 
multiplication number (MN) and, therefore, for the pat- 
tern ([3]), one can write MN = 4. 

The symmetry group G ~ [a*, i] of the bush B[a*, i] is 
determined by two generators: the translation (a^) by 
four lattice spacing (a) and the inversion (j) with re- 
spect to the center of the chain (note that the condition 
A'' mod 4 = must hold for existence of such bush) . 

If we decompose the displacement pattern ^ into the 
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linear normal coordinates [28 



1 



N 

{j = O..N-l), 



27rj \ /27rj 

sm I n + cos n 

N \ N 



n = l..A^ 



(4) 



we get the following form of the bush B[a^, i] in the modal 
space: 



X{t) ^ ^l{t)■^N/2 + J^(t)*3W/4, 



(5) 



where 



vfjv/a = ^{-1, 1 I -1,1 I -1,1 I -1,1 I ...}, (6) 



3^/4 = -^{-1' I'll - 1, -1, 1, 1 I . • •}, (7) 



while and v{t) are time-dependent coefficients in 
front of the normal coordinates m/2 and ^^n/a- 

Thus, only m — 2 normal coordinates from the full 
set (U) contribute to the "configuration vector" X{t) cor- 
responding to the given bush and we will refer to m as 
the hush dimension. 

In 0, we developed a simple crystallographic method 
for obtaining the displacement pattern X[t) for any sub- 
group G of the parent group Gq. Using this method one 
can obtain bushes of different dimensions for an arbitrary 
nonlinear chain. The one- dimensional hushes (m = 1) 
represent symmetry-determined nonlinear normal modes. 
The displacement pattern X{t) corresponding to a given 
NNM depends on only one (time-periodic) function v{t): 



X{t) = v{t) ■ c, 



(8) 



where c is a constant vector, which is formed by the coeffi- 
cients Ci {i = 1..N) from Eq. ([1]), while the function ^{t) 
satisfies a certain differential equation. This so-called 
"governing" equation can be obtained by substitution of 
the ansatz ^ into the dynamical equations of the con- 
sidered chain. 

In some sense, the concept of bushes of vibrational 
modes can be considered as a certain generalization of the 
notion of NNMs by Rosenberg. Indeed, if we substitute 
the ansatz ([5]) into dynamical equations of the chain, we 
obviously get two "governing" equations for the functions 
i'(i) and ii{t), that determines the above-discussed two- 
dimensional bush (note that, in contrast to a NNM, such 
dynamical object describes, in general, a quasiperiodic 
motion). Finally, one can conclude that m-dimensional 
bush is determined by m time-dependent functions for 
which m governing differential equations can be obtained 
from the dynamical equations of the considered mechan- 
ical system. 

Let us emphasize that bushes of modes represent a 
new type of exact excitations in nonlinear systems with 
discrete symmetries and the excitation energy proves to 



be trapped in a given bush for the case of Hamiltonian 
systems. 

It is very important to emphasize that there exist only 
a finite numher of vibrational bushes of any fixed dimen- 
sion in every iV-particle mechanical system. As a conse- 
quence, there is sufficiently small number of NNMs (one- 
dimensional bushes) in the FPU chains (three NNMs for 
the FPU-a model and six — for the FPU-/3 model). 

All possible one-dimensional bushes are explicitly 
listed in our papers 0,0 (see also 

The stability of some NNMs in the FPU chains has 
been studied in [l3 - |22| by numerical and analyti- 
cal methods. Let us comment explicitly on the recent 
papers USUI. 

In [1, [2l|, T. Bountis and coworkers have investigated 
the destabilization thresholds (i?i« and E2u) of two non- 
linear normal modes which they call SPO-1 and SPO-2 
(simple periodic orbits) by numerical methods. The au- 
thors of the above papers try to reveal some relations 
between the destabilization thresholds Ei^ , and the 
origin of the weak chaos in connection with the famous 
Fermi-Pasta-Ulam problem of the energy equipartition 
between different modes. In particular, they conclude 
that the main role in the weak chaos appearance in the 
thermodynamic limit [N — >■ oo) plays SPO-2, because 
E2u ^ Eiu ^ jj, and, therefore, E2u < 

However, there are some other SPOs in the FPU-/3 
chain and one can be interested in their role in the origin 
of the weak chaos in the thermodynamic limit. Some 
comments are appropriate at this point. 

According to Lyapunov , some strictly periodic or- 
bits for nonlinear systems can be obtained from the lin- 
ear normal modes (which are introduced in the harmonic 
approximation) by continuation with respect to a param- 
eter characterizing the strength of nonlinearity. [29] From 
this point of view, there exist N different SPOs for lon- 
gitudinal vibrations of an N particle monoatomic chain. 
However, only few of the modes, constructed in such a 
way, possess an identical time dependence of the displace- 
ments of all the particles. More exactly, only few of the 
Lyapunov modes can be written in the form ([5]) implying 
a separation of time and space variables that is typical 
for the Rosenberg nonlinear normal modes. Indeed, in 
general case, Xi(t) = Vi(t) ■ ct where i>i{t) (i = 1, . . . ,N) 
are different functions of time with identical periods. 

Note that in the present paper we consider only ex- 
tended SPOs, but the same problem there exists for lo- 
calized periodic modes (discrete breathers) and we have 
discussed it in detail in 24] . 

As far as we aware, all periodic solutions in 
monoatomic chains that have been studied up to now 
(see the above cited papers flK^ and references therein) 
belong namely to the class of the Rosenberg nonlinear 
normal modes determined by Eq. Moreover, the spa- 
tial profiles c = {ci, C2, . . . , Cat} of these modes possess 
certain symmetry properties. In particular, every such 
mode can be characterized by a multiplication number 
(MN) determining the enlargement of the primitive cell 
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of the vibrational state in comparison with that of the 
equihbrium state. As was aheady noted, we refer to 
these modes as symmetry-determined NNMs and there 
exist only finite number of such modes (even for the case 

— >■ oo!) for each nonlinear chain [sl [g. [l3|. 

Above considered SPO-1, SPO-2 and the well-known 
TT-mode (zone boundary mode) represent NNMs with 
multiplication numbers 4, 3, and 2, respectively. How- 
ever, among six symmetry-determined NNMs in the 
FPU-/3 chain there exist another three NNMs with 
MN = 3, MN = 4 and MN = 6. 

The stability properties of the second NNM with 
MN = 4 were studied by M. Leo and R.A. Leo in 0. The 
stability of this mode was investigated in the thermody- 
namic limit by both numerical and analytical methods. 

The stability diagrams for all the nonlinear normal 
modes in the FPU-/3 chain, as well as for the FPU-a 
chain, can be found in our paper [6|. With the aid of 
these diagrams, one may reveal many stability proper- 
ties of NNMs for an arbitrary N, in particular, for the 
thermodynamic limit (A^ — >■ oo). Note that these dia- 
grams were obtained numerically. 

In this paper, we present some analytical results for 
the stability properties of all NNMs in the FPU-/3 chain 
in the thermodynamic limit (A^ — > oo). We also com- 
pare our results with those by different authors when it 
is possible. 

In Sec. 2, we consider all the possible symmetry- 
determined nonlinear normal modes in the FPU-/? chain. 
In Sec. 3, the stability diagrams for these NNMs are dis- 
cussed. In Sec. 4, the analytical method for studying 
the stability of NNMs in the thermodynamic limit is pre- 
sented. In Sec. 5, we list our results on the stability 
properties for every NNM. 



II. NONLINEAR NORMAL MODES IN THE 
FPU-^ CHAIN 

As was already mentioned, there exists only finite num- 
ber of symmetry-determined NNMs in any monoatomic 
chain. Every NNM corresponds to a certain subgroup of 
the symmetry group of the chain dynamical equations. 
The difference in the number of nonlinear normal modes 
for the FPU-a chain (three NNMs) and the FPU-/? chain 
(six NNMs) is associated with the fact that the symmetry 
group of the FPU-/? chain dynamical equations is higher 
than that of the FPU-a chain 0, i, d . 

In @, we have investigated the stability of all NNMs 
both in the FPU-a and FPU-/3 chains (for the case /3 > 0) 
by numerical methods. Let us comment on the main idea 
of this investigation. 

Following the standard method of the linear stability 
analysis, we linearize the FPU-/? dynamical equations 
near a given NNM and get the linearized system in the 

form S = J{t) S, where S{t) represents the infinitesimal 
perturbation vector, while J{t) is the Jacobian matrix 
of the original system of nonlinear differential equations. 



Thus, we obtain N linear differential equations with time- 
periodic coefficients depending on the considered NNM. 
Then the Floquet method can be applied for studying 

the stability of the zero solution of the system i5 = J{t) 6. 

However, such straightforward way for the stability 
analyzing becomes practically impossible for N ^ oo. 
In [25l |. the general group-theoretical method has been 
developed for splitting (decomposition) of the original 

system S = J{t) 5 oi N linear differential equations 
into certain subsystems of sufficiently small dimensions 
Nj <^ N . For the FPU-/? chain, these dimensions do not 
exceed three (see below). Then we have applied the Flo- 
quet method for such subsystems of small dimensions. 
Moreover, proceeding in this manner, one can reveal 
those subsets of the vibrational modes, which are respon- 
sible for the loss of stability of the considered NNM. As a 
consequence of this approach, it proves to be possible to 
construct very transparent diagrams, which demonstrate 
explicitly stability properties of each FPU nonlinear nor- 
mal mode 

The FPU-/? model represents a chain of unit masses 
coupled with each other by the appropriate nonlinear 
springs. The dynamical equations describing longitudi- 
nal vibrations of the FPU-/? chain can be written in the 
form 

f{xi+i-Xi)- f{xi-Xi-i), i^l..N, (9) 

where Xi (t) is the displacement of the ith particle from its 
equilibrium state at the instant t, while the force /(Ax) 
depends on the spring deformation Ax as 

f{Ax) ^ Ax + l3{Axf. (10) 

The periodic boundary condition is assumed to hold: 

XN+l{t) = Xi{t), Xo{t)=XN{t). (11) 

Let us mention some results of the paper ,Ql , which are 
necessary for our further discussions. 

Every NNM in the FPU-/? chain can be written as 
follows [see Eq. 

X{t) = iy{t) ■ c, 
where i'{t) satisfies the Duffing equation 



••2 H 
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(12) 



with different values uj and 7 for different NNMs. The 
function vit) describes the time-evolution of a given 
NNM, while the A^-dimensional vector c determines the 
pattern of the displacements of all particles of the chain. 
Below, we list all possible NNMs in the FPU-/? chain. 



1. B[a2,i 



1 



I 1,-1 1 1,-1 1 ...}, 

7 = 16 (/V mod 2 = 0). 
This is a boundary zone mode or 7r-mode. 



;2=4, 



(13) 
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2. B[a^,i]: 

c^^{l,0,-l\l,0,^l\l,0,~l\ ■■■}, .^4. 
^2^3, 7=f (iV mod 3 = 0). ^ ' 

There exist three "dynamical domains" of this 
NNM (see below). 

3. B[a'*, ai]: 

^=7W{0'1'0'-1 I OA,0,~l I 0,1,0,-1 I ...}, 
w2 ^ 2, 7 = 4 {N mod 4 = 0). 

(15) 

There exist two dynamical domains of this NNM. 

4. B[a^,iu]: 

c= ^{1,-2,1 I 1,-2,1 I 1,-2,1 I ...|, , , 
c^2^3, 7=f (TV mod 3 = 0). ^ ' 

There exist three dynamical domains of this NNM. 

5. B[a'^, iu]: 

^2=2, 7 = 8 (A^ mod 4 = 0). 

(17) 

There exist two dynamical domains of this NNM. 

6. B[a'^u, ai]: 

'^^ tIf^O, 1,1, 0,-1,-1 I 0,1,1,0,-1,-1 I ...}, 
cj2 = 1, 7=1 {N mod 6 = 0). 

(18) 

There exist three dynamical domains of this NNM. 

Let us comment on the above listed NNMs in the FPU- 
P chain. 

Every NNM, denoted by the symbol B[G], is charac- 
terized by the corresponding symmetry group G, that 
represents a certain subgroup of the symmetry group 
Go = [a,i,u] of the FPU-/3 dynamical equations (|9I10I) . 
We determine every such group by the set of its genera- 
tors using the following notations: 

a — the translation of the chain by one lattice spacing, 

i — the inversion with respect to the center of the chain, 

u — the operator, that changes signs of the displace- 
ments of all particles without any their transposi- 
tion. 

The symmetry group Go = [a,i,w] of the FPU-/3 dy- 
namical equations is described by three generators (a, 
I, and u). The corresponding transformations a, i and 
u of A'^-dimensional vectors x = {xi,X2, . ■ ■ tXn} do not 
change the dynamical equations (|9ll0p of the FPU-/3 
chain. 

All the above listed groups of NNMs are fully described 
by only two generators, but these generators can be writ- 
ten as some products of the generators a, i, and u of the 



group Go. For example, a^, a^, are translations of the 
chain by two, three and four lattice spacings, respectively. 
The transformation ai means that we must perform the 
inversion of the displacement pattern with respect to the 
chain center and then translate it by one lattice spacing. 
Note that transformations a and i do not commute: 

ia = a~^i, or ia — a^~^i (19) 

[the relation a^^ — a^^^ holds because of the periodic 
boundary condition (fTTj) ]. On the other hand, the trans- 
formation u does commute with both a and i transfor- 
mations: 

ua — au, ui — iu. (20) 

The transformation a^u means that we must change 
signs of all displacements and then translate the displace- 
ment pattern X by two lattice spacings. 

Some simple examples are worth mentioning at this 
point. For the chain with A'^ = 6 particles, we can write 
the following relations: 

a{xi,X2,X3,X4,xz,xe} = {x2,X3,X4,,X5,Xe,xi}, 

i{xi,X2,X3,X4,X5,Xe} = 

{-xe, -X5, ~X4, -xs, ~X2,~'Xi}, 
iu {xi,X2, X3,X4, X5,xe} = {xq, X5,X4, X3, ^2, a;i}, 
a^i {xi,X2,X3, X4, X5,xe} = 

{-Xi, -X3, -X2,-xi,-xe, -X5}, 
a^u {xi, X2,X3, X4, X5, Xe} = 

{-X3, -X4, -X5,-Xe, -Xi,-X2}, 

etc. 

The displacement pattern corresponding to a given 
NNM can be obtained as invariant vector of its sym- 
metry group G C Gq. For example, let us obtain the 
displacement pattern for the NNM with G = [a"' , ai] [see 
Eq. (jl5l) ]. For simplicity, we demonstrate the method for 
obtaining displacement patterns with the case N — 8. 
Let 

X = {xi,X2,X3,X4,X5,Xe,XT,Xs,} 

where Xi {i — 1..8) are arbitrary displacements of eight 
particles of the chain. The vector X must be invariant 
with respect to the action of our two generators and 
ai of the symmetry group of the considered NNM: 

a'^X = X, aiX = X. (21) 

The former equation is reduced to the following form: 

a'^X = {x5,X6,a;7,a;8 | xi,X2,X3,Xi} = 

{xi,X2,X3,Xi \ X^,X(i,XT,Xe,}, 

from which we conclude that 

X^ = Xi, Xq = X2, XT = X3, Xs=X4. 

This displacement pattern is formed by two primitive 
cells whose size four times larger than that of the FPU-/3 
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chain in its equihbrium state. The sets of the displace- 
ments in both cells are identical: 

X {xi,X2,Xz,Xi I Xi,X2,X'i,Xi}. (22) 

Now let us take into account the second equation ([?!]) . 
Acting on the vector (j22p by ai, we obtain 

aiX = a{~XA, -X3, -X2, -Xi \ - x^, -X3, -X2, -xi} = 

= {-a;3, -X2, -Xi, -X4 I - X3, -X2, -Xi,-X4:}. 

Then using the equation aiX = X , we get 

Xi — —X3, X2 — —X2 ~ 0, .T4 = —X4 = 0. 

Thus, the invariant (under the action of the group 
G = [a'*, ai]) vector X depends on only one arbitrary pa- 
rameter, which we denote by x: 

X ^ {x,0,-x,0 \ x,0,-x,0}. (23) 

(Note that this vector being invariant with respect to 
generators of the group G= [a'^,ai] will automatically 
be invariant relative to all its other elements). Then the 
NNM corresponding to the invariant vector (j23l) can be 
written as follows 

XNNAlit) = Mt), 0, ~iy{t),0 I 1/(0, 0, -l^{t), 0} = 

j.(i)-{l, 0,-1,0 I 1,0,-1,0}. (24) 

To find all NNMs, we can try all subgroups of the sym- 
metry group Go = to choose those displacement 
patterns, which depend on only one arbitrary parameter. 

The patterns depending on m arbitrary parameters 
with m > 1 form the m-dimensional bushes of vibra- 
tional modes. Namely in this sense nonlinear normal 
modes may be called one-dimensional bushes. 

In [l,[l^, three different group-theoretical methods for 
constructing the bushes of vibrational modes in arbitrary 
A^-particle nonlinear mechanical systems were developed. 
The most efficient of these methods uses the concept of 
irreducible representations of the symmetry groups. 

Taking into account the above method that was used 
for constructing Eq. we conclude that every NNM 

can be written in the form 

XNNAlit) ^iy{t)-c, (25) 

where cis a certain time-independent vector. Substitut- 
ing ansatz (^5]) into the dynamical equations (PV [TU|) of 
the FPU-/3 chain, with explicit forms of the vectors c 
from Eqs. (fT5HT5|) . one can find that FPU-/3 equations 
are reduced to only one differential equation (governing 
equation of the corresponding NNM) of the form: 

i) + u}'^i' + -f—v^ =0. (26) 

This is the Duffing equation with different values uj and 
7 for different NNMs which are hsted in Eqs. PMT5)) . 



Above, we have mentioned the existence of so-called 
"dynamical domains" of all nonlinear normal modes pre- 
sented in Eqs. (jl3ffl8l) . Let us comment on this notion 
borrowed from the theory of phase transitions. 

We have already emphasized that a certain symme- 
try group G corresponds to every NNM. This group is a 
subgroup of the symmetry group of the considered me- 
chanical system in its equilibrium state (G C Go). If we 
act on the vector X{t) corresponding to a given NNM by 
operator g £ Gq, that does not belong to subgroup G, 

we get the equivalent configuration vector X{t) — gX{t). 

The equivalent vector X{t) corresponds to a new NNM, 
which is described by the same dynamical equations as 
that of the NNM associated with the vector X{t). For 
example, three dynamical domains are associated with 
the NNM from Eq. ([Ti)) : 

B[a3,z]:c= -£={1,0,-1 | 1,0,-1 | 1,0,-1 | . . .}, 
VoA' 

(27) 

B[a^at] :c= -£={0,1,-1 | 0,1,-1 | 0,1,-1 | . . .}, 
VoA' 

(28) 

B[a^a2^] :c= —£={1,-1,0 | 1,-1,0 | 1,-1,0 | . . .}. 

(29) 

All the displacement patterns (PTH^^ differ from each 
other by a cyclic transposition of the displacements inside 
each primitive cell of the chain equilibrium state. Let us 
note that the symmetry groups Gj [j = 1,2,3) of NNMs 
from Eqs. (|27H29p prove to be conjugate subgroups in the 
parent group Go, for example, G2 = g~^Gig {g £ Go). 

Since the above-discussed "domains" possess equiva- 
lent dynamical properties, we study below the stability 
of only one copy of the full set of dynamical domains for 
every NNM in the FPU-/? chain. 

All symmetry-determined NNMs that can exist in the 
FPU-/3 chain with an appropriate number of particles are 
listed in Table HI 

III. STABILITY DIAGRAMS FOR NNMS IN 
THE FPU-/3 CHAIN 

The detailed numerical analysis of stability of NNMs 
in the FPU chains can be found in [6|, where results 
have been presented in the form of certain "stability di- 
agrams" . Let us consider the structure of such diagrams 
and the method of their obtaining. 

Every NNM represents a periodic dynamical regime 
and the standard Floquet method can be applied for in- 
vestigation of its stability. According to this method, we 
have to linearize the nonlinear FPU-/3 dynamical equa- 
tions (|9ll0p in the vicinity of a given NNM XNNAiit)- 
To this end, let us introduce an infinitesimal vector 

6 = {Si,62,S3,...,Sn}, (30) 
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TABLE I: Nonlinear normal modes in the FPU-/? chain 



NNM 


Displacement 


Modal space 


Governin 


B equation 




pattern 


represent at ion 








B[a^,i] 






i> 




if ^.3=0 

iV 


B[a^i] 


|i,-i,o| 


I^^(*JV/3 - *2JV/3) 


i/ 




2N ^ 


B[a^,m] 


|1,1,-21 


I^75(1'iV/3 + *2iV/3) 


V 




1^.^ = 


B[a^, ai] 


11,0,-1,01 


«^75(*iV/4 - *3JV/4) 


i) 


+ 2p + 




B[a'^,iu] 


11,-1,-1,11 




i) 


+ 2v + 


f^^^o 


B[a^u, ai] 


11,1,0,-1,-1,01 


I^73(*iV/6-*5iV/6) 


V 







that determines a certain perturbation of the NNM. Let- 
ting 

x{t) = XNNMit)+m, (31) 

we substitute X{t) into nonlinear equations (j9ll0p and 
omit all nonlinear terms in 5j{t) {j — 1..N). As a result of 
this procedure, we get the linearized (variational) system 

5 = J{t) 6, (32) 

where J{t) is the corresponding Jacobian matrix (see de- 
tails in [fjj and especially in [25]). 

Equation ((32)) represents an TV x system of linear 
differential equations with time-dependent coefficients, 
which, in turn, are determined by the periodic func- 



tion ^{t) describing the time-evolution of the considered 
NNM. 

In [2^ , we have developed a group-theoretical method 

for splitting the linearized system S = J(t) 6 into a num- 
ber of independent subsystems of small dimensions (for 
all NNMs in the FPU-/? chain these dimensions do not 
exceed 3). In the explicit form, this splitting can be 
found for all NNMs in the FPU-/3 chain in Table 8 of 
the paper Let us comment on that table using as an 
example the nonlinear normal mode B[a'*,m] for /3 > 0. 
The "rabbit ears" stability diagram, depicted in Fig. [TJ 
corresponds to this mode. 

In this case, the linearized system can be splitted into 
(^ — 2) two-dimensional systems of the following form 



-4sinM^ 



5,,+4sinM ^ 



6/3 2/ 



12/3 



^^^(i)sin — iU, 



N 

-if.^(t)sin(^)5, 



2Trj 



N 
N 



(33) 



Here f = f - J, and j = 1, . . . , (f - l) . 

Thus, the jth component of the infinitesimal vector 6 
(after the appropriate orthogonal transformation) turns 
out to be coupled only with its j' component, where 
j' = Y ~ j ■ Therefore, the investigation of stability of 
the nonlinear normal mode B[a^,m] reduces to studying 
stability of zero solution of Eqs. (j33p with different val- 
ues of the indices j and j' ~ ^ — j ■ It can be seen from 
Eqs. ([55]) that the pairs of modes whose indices are situ- 
ated symmetrically with respect to the index jo — ^ oi 
the considered NNM are coupled. 

Note that Eqs. ([55]) represent the system of two cou- 
pled differential equations with time-periodic coefficients 
determined by the solution v{t) of the Duffing equa- 
tion (pS)) . It is well known that equations of such a type 
possess domains of stable and unstable motion, [s^ 

For the system (j33p . we can construct the stability di- 
agram in the plane (j — A), where A is the amplitude of 



the NNM [i.e. the amplitude of the function v{t) from 
Eq. ([26])], while j is its index. Such diagrams provide us 
with the most valuable information for the case of the 
"thermodynamic limit" when the number of lattice cells 
tends to infinity {N -^00). 

When A^ — >■ 00, the number of modes also tends to in- 
finity and it is convenient to introduce the wave number 
k {{) < k < 2tt) instead of the mode index j: 

J^l-N. (34) 

In Eq. (HI there is a normalizing coefficient 1 / \/N origi- 
nating from the conventional definition of the mode norm 
I'fjP = 1. As a consequence, displacements of all par- 
ticles corresponding to a given normal mode with fixed 
amplitude tend to zero with A^ — > 00, while the coeffi- 
cient 1/A^ appears in Eqs. (l26l) . ([33]) etc. in combination 
with the coefficient /3. 
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Hereafter, we use a more convenient for our purposes 
normalization of normal modes (normalization to the 
"volume" of the system), namely, we assume IvPjP — N. 
This permits us to exclude the coefficient 1/N coupled 
to /3. Moreover, using a trivial coordinate rescaling, one 
can set \f3\ = 1. 

As a result, Eqs. (1551) are transformed to the following 
form: 

4+4sin2(|)[l±6^2(^)]j^ ^ Tl2j^'(t)sin(fc)4', 
4' +4cos2(|)[l±6i/2(t)]4, = zfl2„^{t)sm{k)Sk. 

(35) 

Here the upper and lower signs correspond to /? > and 
/3 < 0, respectively. 

Following this idea, we construct our stability diagrams 
in the (fc — A) plane, where A is the amplitude of a nor- 
mal mode subjected to the above mentioned normalizing 
condition jvl'jp = N. For any fixed N, the permissible 
values of the wave number k represent the equidistant set 
inside the interval [0, 27r]. When N oo, these values of 
k form a dense set on the above interval, depicted on the 
horizontal axis of diagrams in Fig. [1] 

Rich information about stability properties of NNMs 
can be obtained from diagrams of such a type. Let us 
consider this idea in more detail using as an example the 
stability diagram of the NNM B[a*, iu] depicted in Fig.[T]3 
for /3 > 0. 

In this diagram, we denote the permissible values of the 
wave number fc for iV = 12 by the dotted vertical lines 
{k = 7, 27, 37, ... , 127, where 7 = 27r/12). The black 
color corresponds to the regions of unstable motion in 
the plane {k — A). From Fig. [T^ one can see that for the 
FPU-/? chain with = 12 particles the NNM B[a^,iu] 
with c = {1, -1, -1, 1 I 1,-1,-1, 1 I 1, -1, -1, 1} turns 
out to be stable for every amplitude A up to 2 (actu- 
ally, at least, up to A = 20). Indeed, the dotted vertical 
lines do not cross the black (unstable) regions in the form 
of the rabbit ears. 

Obviously, the mode B[a'',iu] is also stable for the 
chains with = 4 and A^ = 8 particles since the verti- 
cal dotted lines, similar to those depicted in Fig. [T^, are 
more distant from each other than for A^ = 12, and the 
regions of instability situate fully between the neighbor- 
ing lines. 

In the case A^ = 16, there exist such dotted verti- 
cal lines that are tangents to the rabbit-ears unsta- 
ble regions, while for A^ > 16 (note that the relation 
A^ mod 4 = must hold!) these lines begin to cross the 
unstable regions. Therefore, we can conclude that the 
considered nonlinear normal mode in the FPU-/3 chains 
with A^ > 16 particles becomes unstable for the vibra- 
tional amplitudes, which fall into the black regions in 
Fig. [T]3. The case A" = 16 represents the "boundary" 
case between stable and unstable behavior of the B[a'*, iu] 
NNM for appropriate amplitudes. 

From Fig.[T^, it can be seen that the critical amplitude 
Ac of the nonlinear normal mode B[a'*, iu], for which this 
mode loses its stability, decreases with increasing A^ and 



we can evaluate numerically the corresponding scaling 
law of the function Ac{N) in the thermodynamic limit 
N 00. 

Let us consider this question in more detail. The 
instability regions touch the fc-axis at the points ^ 
(j = ^) and ^ (j = ^) , which correspond to the NNM 
B[a'*,m] and its dynamical domain B[a'*,a^m], respec- 
tively. Namely near these points the considered NNMs 
loses their stability in the thermodynamic limit N 00. 

Let us focus on the point (fc = which corresponds 
to the index jo = -j of the mode B[a^, iu]. AU the neigh- 
boring fc-points correspond to the "sleeping" modes, i.e. 
modes with zero amplitudes for the stable behavior of 
the mode B[a^,iu]. If we increase the amplitude A of 
the NNM B[a'*,m] from zero, the sleeping modes, which 
become firstly excited, are those with indices j' = ^ ± 1, 
i.e. the modes closest in number to the index jo = "j '-'■^ 
the NNM B[a^,iu]. 

We denote by Ac{N) the critical value of the amplitude 
A of the considered nonlinear normal mode for which the 
loss of its stability takes place, i.e. Ac{N) is the threshold 
of the B[a'', iu] NNM stability for a given A^. Note that 

we can speak about the threshold value Ai^^ for the exci- 
tation of the sleeping j-mode as a result of its interaction 
with the B[a'*,ju] NNM whose index is equal to jo = "j. 
Then from Fig. [TJs, we obtain 

Ac{N) = min^^') = A^c^'^^\ (36) 
3 

For the index j situated near jo — ^ ~ f ) > '^^'^ '-^^^ 
see from Fig. [TJ?, that 

= Afc -tana, (37) 

where Ak = ^ is the minimal distance between the 
neighboring permitted points on the fc-axis and tana is 
determined by the tangent to the black unstable region 
near the point fco = f • 

The scaling law of the function Ac{N) for A^ — )■ 00 can 
be deduced from Eqs. ([57)) : 

Ac(A) = ^tana {N ^ 00). (38) 

Here tana '-^ 1, as it might be evaluated from Fig. [TJj by 
the straightforward manner. 

The above discussed threshold function Ac{N) for 
A^ — > 00 can be calculated analytically. For the nonlinear 
normal mode B[a^, iu], such calculation has been recently 
performed in [21 . The similar analytical results for the n- 
mode (jo = ^) were earlier obtained in [H. [l^ [l6l - l20j . 

In the present paper, we present the analytical estima- 
tions of stability thresholds for N ^ 00 for all possible in 
the FPU-/3 (/3 > and /3 < 0) nonlinear normal modes 
(see Sec. 5). 

Let us return to the stability diagram in Fig. [T]; for 
the NNM B[a'*,m] for (3 > 0. It is very interesting that 
for the sufficiently large amplitudes (A > Ac ~ 0.346) the 
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f) 








27r 



FIG. 1: Stability diagrams for NNMs in the FPU-/3 chain (left ,8 > 0, right /3 < 0): a) B[a^,i]; b) B[a^,i]; c) B[a*,ai]; 
d) B[a^,iu]; e) B[a'',w]; f) B[a^ u, ai]. 
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nonlinear normal mode B[a^,iu] again becomes stable 
even in the thermodynamic limit N ^ oo\ Results of 
such a kind cannot be obtained by the analytical meth- 
ods applied in 2] and in the present paper. However, 
the stability of NNMs for nonlinear chains in the case 
of large amplitudes can be analyzed by the numerical 
method used in p'l which provide us the stability dia- 
grams depicted in Fig. [T] 

In the recent paper [l^l, the existence of the "second 
stability threshold" Ec for the NNM B[a^,m] was also 
revealed. Beyond Ec, the stability of the mode B[a'*,m] 
(in [l^l, it is called 7r/2-mode) is restored. This result 
was obtained by numeric methods. In term of the partial 
energy per one particle of the FPU-/3 chain, this threshold 
turns out to be equal to Ec = 0.14715. On the other 
hand, we can find an approximate value for Ec directly 
from diagram depicted in Fig.[T]3: Ec ~ 0.15. This value 
is in good agreement with that from 22]. 

Moreover, such a "surprising behaviour" (as has been 
written in [22|) one can reveal directly from our sta- 
bility diagrams presented in for NNMs B [a^ , i] and 
B[a'^,iu]. We reproduce these diagrams in Fig. [TJd and where 



[TJi of the present paper, and one can found from them 
the approximate values of the "second stability thresh- 
olds" for the above mentioned nonlinear normal modes: 
Ec{B[a^,i]) « 3.70 and Ec{B[a^,iu]) « 24.8. 

In the right column of Fig. [U we present new results 
on the stability of all possible NNMs in the FPU-/3 chain 
with /? < 0. In contrast to the case /3 > 0, the chain 
with (3 < occurs to be unstable itself for the mode 
amplitudes which exceed a certain critical value (this 
case corresponds to the region depicted by gray color 
in Fig. H]). Note that beyond E = 0.092 for the NNM 
B[a^,i] and E = 0.14 for the NNM B[a^,iu], all sleeping 
modes turn out to be excited simultaneously. Such a re- 
markable behaviour we have previously revealed for the 
stability properties of the FPU-a chain (see [1]). 



IV. ANALYTICAL METHOD FOR STABILITY 
ANALYSIS OF NONLINEAR NORMAL MODES 
IN THE THERMODYNAMIC LIMIT 

For studying stability of NNMs in the thermodynamic 
limit {N — > oo), we use a method which is similar, 
in some sense, to that developed in [2]. However, our 
method is more general and, as a consequence, it can be 
applied to all nonlinear normal modes in the FPU-/3 chain 
for the case N — ?> oo. Below, this method is described in 
detail. 

The standard linear stability analysis of a given NNM 
leads us to Eq. (15^ representing N x N system of dif- 
ferential equations with periodic coefficients. Obviously, 
such straightforward way proves to be especially compli- 
cated when iV — >• oo. 

On the other hand, we can apply the general group- 
theoretical method [23 for splitting the system ([32]) into 
some independent subsystems Lj of small dimensionali- 



ties nj due to the symmetry properties of the considered 
NNM. In it has been shown that for all NNMs in the 
FPU-/3 chain these dimensionalities {nj) do not exceed 
3. This fact is extremely useful for stability studying and 
we have already used it in [6] for numerical construction 
of the stability diagrams for NNMs in the FPU-/? model 
with 13 >0. 

At this point it is appropriate to dwell on the physi- 
cal cause of the stability loss of nonlinear normal modes. 
This loss of stability can be treated in terms of "para- 
metric interactions" between a given excited NNM and 
other modes with zero amplitudes ("sleeping" modes) 
(see [IH). 

In the simplest case, such interaction can be described 
by the Mathieu equation. Indeed, studying stability of 
the TT-mode in the FPU-/3 chain, we deal with splitting 
of the A^-dimensional variational system into individual 
scalar equations 



1 



12;3 



0, 



(39) 



Vj{t) is a sleeping mode, while the solu- 



tion i'{t) of the Duffing equation ([26)) describes the time- 
evolution of the TT-mode. 

In some approximation, we can replace the exact func- 
tion i^(t) by its first Fourier harmonic. After substitu- 
tion of this approximate function into Eq. (I39p . the latter 
can be transformed to the standard form of the Mathieu 
equation, which possesses an infinite set of stable and 
unstable regions. Depending on the value of the con- 
sidered NNM amplitude [the function ^{t)] one can get 
into a region of unstable motion and then infinitesimal 
solution of the Mathieu equation begins exponentially in- 
crease in amplitude over time evolution. This means that 
the original NNM loses its stability and transforms into 
a two-dimensional bush of vibrational modes. 

The above phenomenon represents a parametric res- 
onance: the parameter in square brackets in Eq. (j39p 
changes periodically that leads, under certain conditions, 
to excitation of the sleeping degree of freedom described 
by the function Vj{t). 

Let us note that different sleeping modes i'j{t) turn 
out to be unequal in relation to excitation by the given 
NNM v{t). Indeed, in Fig. [TJi it can be seen that Vj{t) 
with wavenumbers k closer to fcrcs = 7r/2 prove to be 
excited earlier, i.e. for smaller amplitudes A of the tt- 
mode. The values k^es for different NNMs in the limit 
N ^ CO can be revealed in Fig. [TJ they correspond to 
the points that contact with the fc-axis. We need to study 
vicinities of the points (fc = fcres, A = Q) for obtaining 
scaling of the stability thresholds in the thermodynamic 
limit (A'' — >■ oo). 

In general case, parametric interaction with the origi- 
nal (active) NNM leads to exciting not only one sleeping 
mode, but a certain set of such modes (for the FPU-/3 this 
set can consist of 1, 2 or 3 sleeping modes). For example, 
the variational system for NNM B[a^,m] can be decom- 
posed into two-dimensional subsystems [see Eq. ((33)) 
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and the discussion below this equation]. Each of these 
subsystems contains two sleeping modes Vj (t) and Vji (t) , 
where j' — y ^ ^ terms of wavenumber (j34p . these 
modes are denoted by i'k{t) and v:^-kit)]. This means 
that when NNM i'{t) loses its stability the both sleeping 
modes i^fc(i) and v^L^kit) are excited simultaneously. 

Similarly, the variational system for the NNMs B[a'^, i] 
and B[a'^,m] can be decomposed into three-dimensional 
independent subsystems ■ This means that these 
nonlinear normal modes lose their stability simultane- 
ously relative to the three sleeping modes with wavenum- 
bers k,k + ^-,k + ^. 

Since there are no analytical results on parametric res- 
onance in such cases, we have to use the general Floquet 
stability analysis to obtain thresholds of the stability loss 
of the corresponding NNMs. 

In terms of the Floquet method, the stability loss of a 
given NNM can be interpreted in the following manner. 
When we vary the amplitude A of NNM in the stabil- 
ity region, Floquet multipliers move on the unit circle. 
Some of them can leave this circle after colliding with 
each other, and this phenomenon means that the consid- 
ered NNM loses stability. On the other hand, Floquet 
multipliers represent eigenvalues of the monodromy ma- 
trix X(7r), and they can be found as the roots \j of its 
characteristic polynomial. 

A bifurcation from stability to instability of the con- 
sidered NNM takes place when some roots of this polyno- 
mial coincide with each other. This fact can be revealed 
by vanishing of the discriminant D of the monodromy 
matrix characteristic polynomial. Indeed, in general case. 



D = WiK - \j 



(40) 



and coincidence of every pair of eigenvalues (Ai, Xj) leads 
to the relation D = 0. 

Now let us consider some technical details of our 
method for studying stability of NNMs in the thermo- 
dynamic limit. 

The explicit form of the independent subsystems Lj 
for all FPU-/3 NNMs can be found in [6] . In the present 
paper, we write Lj for each NNM in the following matrix- 
vector form: 



+ bv^{t)M fi = {) 



(41) 



Here uP' and M are Uj x Uj constant matrices depending 
on the harmonic frequencies of the variables coupled by 
Eq. gl]). 

In Table ini we present the matrices and M for each 
of six NNMs permissible in the FPU-/? chain in terms of 
the frequencies = sin(|), W2 — sin(|-|-|), 0)3 = 
sin(| + f),^4 = sin(| + f). 

The time-depending function v{t) entering Eq. (j4ip is 
a periodic solution to the governing Eq. ([T2|). For ev- 
ery NNM in the FPU-/? chain, this governing equation 



represents Duffing equation 

i) + av + hv'^ ^ 0, 



(42) 



where a is the squared frequency of the harmonic approx- 
imation, while & = • 7 is a nonlinearity coefficient. 

Eq. is called the hard {soft) Duffing equation if 
6 > (fo < 0). For initial conditions z/(0) = A, i>(0) = 0, 
the solution of the hard Duffing equation can be written 
in the form 



Here 



v{t) = Acn(m,K^). 



= a/(l - 2^2), 



(43) 



(44) 



while modulus k, of the Jacobi elliptic cosine is deter- 
mined by the relation 



2k^ ^bA^/{a + hA^). 



(45) 



The solution (|43p represents periodic function with the 
period 



T = 4:K{K)/n, 



(46) 



where K{k) is the complete elliptic integral of the ffist 
kind. 

For the same initial conditions, the solution of the soft 
Duffing equation {b < 0) can be written in the form 



with 



v{t) = Asii{nt + K{k),k^), 



=0/(1 + ^2), 



= \b\A^/{2a-\b\A^), 



(47) 



(48) 



(49) 



T^4:K{K)/n. (50) 
It is convenient to introduce the time scaling 



_ 27r _ irft 



t 



(51) 



which transforms Eq. (I4ip into the equation with tt- 
periodic coefficients. As a result of this scaling the 
form of Eq. (HT|) does not change, but the constant a 
and 6, entering this equation, must be multiplied by 
4:K^{K)/{n^Q'^). We do not change notations in Eq. (|4Tt . 
however, we imply below that the above transformations 
are already fulfilled. 

Our further stability analysis of NNMs reduces to in- 
vestigating the stability of the zero solution of Eq. (j4ip . 
The analysis consists of the following steps: 

Step 1. Simplification of Eq. (|4T]) in the thermody- 
namic limit (A^ — ^ 00). In this case b ^ 1/N and we can 
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TABLE II: Input data for stability analysis of nonlinear normal modes in the FPU-/? chain 



NNM 



Energy vs. NNM's amplitude 



M 



B[a^^] 

B[a^,m] 

B[a*, ai\ 
B[a*, iu\ 

Bfa'^ii, ai\ 



2A^ + f 



3 42 I 27^J 44 
2^ + SlV ^ 



3 42 , 27£ .4 
2^ + 8JV^ 



I 3^44 



2 



2 



CJ2 



a;3 



CJ2 



V 



^3 



aj4 



L02 




UJIUJ4, 



-2 



LU3 



-LOlU]2 



2(1)2 — £l'2<i3 



decompose the coefficients of Eq. (pi]) into power series 
with respect to the small dimensionless parameter 



5A2 



(52) 



Step 2. Search for the general solution of the approx- 
imate equation that was obtained as a result of Step 1. 

Step 3. Construction of the monodromy matrix with 
the aid of the above solution. 

Step 4. Construction of the characteristic polynomial 
of the monodromy matrix. 

Step 5. Analyzing discriminant of the characteristic 
polynomial in the limit N ^ 00. 

Let us consider these steps in turn. 



sn(u, K^) 



2tt 



-Y- 



(2n — 1)ttu 
sm ^ { , (54) 



kK{k) 1 - g2«-i 2K{k 



where 



(2n- 1)!! 



1 2 



2"n! 



(55) 



is the complete elliptic integral of the first kind, while 

K{k) 



q = q{K) = cxp — TT 



Step 1 

The function //(<) in the form for the case f3 > 
and in the form (HT)) for /? < must be substituted into 
Eq. (j41l) taking into account that modulus k goes to zero 
when N 00. To simplify i^(t), we use the following 
formulas from the theory of elliptic functions [1^ : 



cn(u, K^) 



2tt 



00 



kK{k) ^l + q 



2n~l 



COS 



{2n — 1)ttu 
2K{k) ' 



(53) 



(k' = Vl — is the complimentary modulus of the el- 
liptic functions). Note that the modulus n depends on 
hA^ in a different manner for the cases /3 > and /3 < 
[see Eqs. P5|) and (H^ . respectively]. 

Now, we have to decompose the left-hand-side of 
Eq. (|^T|) into the power series with respect to the small 
parameter e = hA^ /a. This very cumbersome decompo- 
sition has been performed with the aid of the Maple^" 
mathematical package. The corresponding result can be 
written as follows: 





3,2 












- - cos 2r I M 








2 ) . 



' 75 . 
.128" ^ 

243 . 

CJ" 

512 
O(e^)=0, 



13 3 „ 1 \ 

-32-r°^'"+32^°^4V''' 
87 597 „ 3 , 
256 + 2048 ^°^ 2" - 64 '^"^ 4" 



(56) 



2048 



cos6t M 
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step 2 



formal series 



Eq. ([55)) represents a system of differential equations 
with time-periodic coefficients and to construct the corre- 
sponding monodromy matrix we must obtain its solution 
for t = TT. On the other hand, for small time intervals, 
the solution of Eq. (j56p can be found by a simple per- 
turbation theory. To that end we decompose fl{t) into a 



substitute it into Eq. (j56|) and equate to zero the terms 
with every fixed power of the small parameter s. As a 
result, we get the following set of differential equations 



flo+uj Mo = 0, 



cos 2r j M 
cos 2t I M 



Ail 



(58a) 
(58b) 

(58c) 



75 
128" 



13 3 „ 1 , , 

cos 2r H cos 4t M 

32 8 32 ' 



Because of the diagonal form of the matrix (see 
Table HI]) . these equations determine certain sets of har- 
monic oscillators with different time-periodic external 
forces. Each of these oscillators is described by equation 

X + uj'^x — /(r). 

The general solution to this equation, obtained by the 
method of variation of arbitrary constants, can be written 
in the form 

T 

x{t) — cis'muT + C2COSUJT -\ / f(t)s'm[uj(T — t)]dt, 

w J 



(59) 

where f{t), in our case, represents a sum of time-periodic 
functions whose frequencies are incommensurable. In- 
deed, for the most NNMs from Table [TTl uj'^ are matri- 
ces with different diagonal elements and, therefore, the 
components of the vector 11q{t) from (|58ap vibrate with 
different frequencies. Substituting /io(r) into (|58bl) leads 
to mixing its time-depended components because of mul- 
tiplying by the matrix M , and such a mixing produces 
more and more complicated terms in r.h.s. of Eqs. (j58p 
when we take into account higher orders in the decom- 
position ([55]) 

Step 3 

The usual way to study stability of a given periodic dy- 
namical regime is the Floquet method. In this method, 
we linearize nonlinear equations of motion in the vicin- 
ity of the periodic solution and calculate the monodromy 



matrix X(T) by integrating 2n times the linearized equa- 
tions with time-periodic coefficients over one period T 
using specific initial conditions [n is the number of equa- 
tions in (155)) ]. These conditions are determined by the 
successive columns of 2n x 2n identity matrix. 

Solving Eqs. (|58p in step-by-step manner, we can con- 
struct the approximate analytical solution to Eq. (j56p up 
to a fixed order of the small parameter e. With the aid 
of this solution, we are able to obtain the monodromy 
matrix X{tt) for Eq. ([55)) . where tt is the period of its 
coefficients. 

The stability of the considered periodic solution is de- 
termined by Floquet multipliers representing eigenvalues 
of the monodromy matrix. If all these multipliers lie on 
the unit circle, the solution is linear stable. In other case, 
the solution linear unstable. 



Step 4 

We obtain eigenvalues of the monodromic matrix X{Tr) 
as the roots of its characteristic polynomial. 

Let us remind that according to the Newton formulas, 
the coefficients of the characteristic polynomial 

det{A-\E) = (-l)^(A^+piA^-i + ---+pAr_iA-|-pw) 
of any N x N matrix A can be expressed via the sums 

JV 

s,-^Af, fc-l,...,iV (60) 

i=i 

with the aid of the recurrence relation 

Pk = -(Sfc +PlSk-l +P2Sk-2 H \-pk-lSl)/k. 
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Thus, we have 



- (S2 +PlSl)/2, 

- (S3 +P1S2 +P2S1), 



Pi = 

P2 = 
P3 = 



On the other hand, all sums Sk in Eq. ([60]) can be found 
directly by means of traces of the matrix A: 



Sk 



Tr(i'=) 



fc = 1, 



.,N. 



It is well known, that in the case of any Hamiltonian 
system with n degree of freedom Floquet multipliers Xi 
form pairs Xi, (i = 1, . . . , N) and, as a consequence, 
the characteristic polynomial /(A) of the monodromic 
matrix X{tt) proves to be palindromic: 



\2n-2 



with the following coefficients 
ai ~ Tr X{it), 



. + a2A^-aiA+l = 0. 

(61) 



1 



Tt^ X{7r) - TrX^in) 



Now we have to obtain formulas for discriminants Di , 
D2 and D3 for the corresponding palindromic charac- 
teristic polynomials of n = 1,2,3 degrees via traces of 
monodromy matrices. 

From the above Newton formulas applied to the poly- 
nomial (j61l) . we obtain the following formulas for its co- 
efficients: 



ai = Tr X(7:), 



a2 = 



0-3 



Tr^X(7r) - TiX^i?:) 



(62) 



Tr^X(7r) - 3TrX(7r)TrX2(7r) +2TrX3(7r) 



On the other hand, one can express discriminant Dn 
explicitly via these coefficients of the characteristic poly- 
nomial. With the aid of Maple'^", we finally find 



sides of the black region in Fig. [T^, which are not excited 
by parametric interactions with the 7r-mode. This fact 
was revealed analytically in [l^ with the aid of Rotating 
Wave Approximation (RWA). In Ts'], in the framework 
of the same approximation the following relation between 
the amplitude threshold value Ac and the wavenumber k 
(i.e. the boundary curve of the black region in Fig. 
was obtained: 



Ar 



sm 



k — 7T I 



9 cos2 



This analytical formula is in good agreement with the 
numerical results. 

Let us now consider the stability threshold of the tt- 
mode in the thermodynamic limit TV — )■ 00 using the 
above discussed method. We have to consider the vicinity 
of the point (f ,0) on the (fc — A) plane (/cj-os = ^f^)- 
The one-dimensional constant matrices of the decoupled 
variational system can be found in Table |lTJ 

w?, = sin2^, Mfc=3sin2^. 

The monodromy matrix also depends on the wavenum- 
ber k and we can decompose its trace TrAr(7r) into the 
Taylor series in two small parameters Ak = k — k-^cs and 
e = hA^ /a (in our case, a — 4, b = 16/3/7V and, therefore, 
e = A/SA^ /N). This decomposition read 



TrX{7:) = ( -2+—TT'Ak 



1536 



32 512 



27 
256' 



n^Ak^e 



2^2 



In the considered case, n = 1 and the corresponding 
discriminant is 

Di = (ai +2)(ai -2), 

where ai ~ TtX{tt). 

The condition Di ~ leads to the equations 



Di =(ai+2)(ai-2), 

D2 =(02 + 2ai + 2){a2 - 2ai + 2)(8 + aj ~ 403)^ 
L>3 =(03 + 2a2 + 2ai + 2)(a3 - 2a2 + 2ai - 2)- 



{9al 



54aia3 — 2703 — A2a\a2 + 18010203 



4(02 - 3)^ + 8af + a\a\ - 4a\a-if 



(63) 
(64) 

(65) 



where Oi, 02, 03 are given by Eqs. ([5^ . 

Now we present some illustrations of the above dis- 
cussed technique. 

Example I. Nonlinear normal mode B[a^,i] (7r-mode): 
\x,-x\. 

Firstly, let us discuss the case /3 > 0. One can see that 
there exist modes, corresponding to the left and right 



oi + 2 = or ai - 2 = 0. (66) 
From the first of these equation, we find 



32 512 



and then we obtain 



e = ^Afc^ ^ —Ak^ ^ O(Ak^). 
6 144 
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Substitution of e 
following result 

TT 

Ar = 



§Al and Afc 



^ leads us to the 



+ 0(iV-9/2) 



The corresponding energy per one particle is 



N 



1 

TV 



2 



0(iV-^). 



The second equation leads to a contradiction with 
the condition of smallness of the parameter e and, there- 
fore, this equation doesn't produce instability of the tt- 
mode in the limit N ^ oo. 

Note that the analytical dependence Ec/N — 
TT^ / (3/?iy) was revealed in [l^, and later was recovered 



in [l|,|Tl 



In the case /3 < 0, the stability properties of the tt- 
mode are utterly different. Indeed, one can see in Fig.[T^ 
(right column) that this mode turns out to be stable up 
to the finite value of its amplitude Ac- Using a numer- 
ical method, we found that in the thermodynamic limit 
TV ^ oo Ac = 0.393 {Ec = 0.213). 

Example 2. Nonlinear normal mode B[a^,m]: 

I X J tC J X ^ X I . 

The variational system for this NNM is decoupled into 
2x2 subsystems with time-periodic coefficients. As a 
result, for studying the stability loss of the mode B[a^, iu] 
we have to vanish the discriminant D2 from (j64p . 

The first factor of D2 near the resonance wavenumber 
kj-cs = 7r/2 is equal to 



a2 + 2ai + 2 



128 



16 



(67) 



Being positive, this factor cannot lead to the condition 
D2 = 0. 

The second factor of D2 from Eq. ([M)) reads 



a2-2ai + 2 = 16 + 0(e2,Afc2) 

also doesn't vanish in the thermodynamic limit iV 
Only the last factor of the discriminant D2 



64 



1 

16' 



can lead to fulfilment of the condition D2 = 0- 
yields 

e = ±-Ak^, 
3 



(68) 
00. 

(69) 
This 



and, therefore. 



Ar 



Er 



3|/3|7V2' 



(70) 



Note that this NNM has been investigated quite 
enough in Above, we simply reproduced the main 
result of this paper by our method. 



Example 3. Nonlinear normal modes B[a^, i]: |x, 0, —x\ 
and B[a"^, iw]: \x^—2x,x\. 

The variational systems for these NNMs can be decou- 
pled into 3x3 independent subsystems whose matrices 
are presented in Table [Til 

Now one has to vanish the discriminant from 
Eq. (l65|) . With the aid of our method, we get the fol- 
lowing results. 

For both modes B[a^, i] and B[a^, iu], we have obtained 
identical scaling for the case N ^ oo: 



Ar.^ 



2tt 



Er. = 



2ti' 
3/3iV2 



(71) 



Some numerical results on stability of the NNM B[a^, i] 
have been found in Q , but we don't know any results on 
stability properties of the NNM B[a'^,m]. 

In conclusion, let us consider another scenario of the 
stability loss of NNMs. Indeed, up to this point, we have 
discussed only the loss of stability associated with para- 
metric interactions of a given NNM with other (linear) 
normal modes of the FPU-/? chain. 

Some NNMs in the FPU-/? chain, when e = 
bA^/a —T' 0, transform not into one linear normal mode 
(LNM), but into a certain superposition of such modes. 
For example, NNM B[a^, ai] 



X{t) = i.{t) 



/2N 



{0,1,0,-1|0,1,0,-1|...} (72) 



transforms, in the case e -> 0, into the linear combination 
ii^N/i{t) + 6*3W/4(i), 6 = -6 = ^ (73) 
of two linear normal modes 



*A'/4(^) = Cjv/4 ■ COs(wi), 

^3N/i{t) = C37V/4 • cos(a;i). 



(74) 



where 



Cn/4 
C3N/4 



^{1,-1,-1,1|1,-1,-1,1|...}, 

^{-1,-1,1,1| -1,-1,1, 1|...}, (75) 
2. 



One can also say that NNM B[a^, ai] is the result of the 
continuation of the superposition (j73p with respect to the 
nonlinearity parameter e of the considered FPU-/? chain. 
We have to emphasize that the continuation of an arbi- 
trary linear combination of two above discussed LNMs 
doesn't represent an exact solution to nonlinear dynami- 
cal equations of the FPU-/? chain, while the superposition 
with ^1 = —^2 = produces an exact solution. 

Then, the following question arises: If we will increase 
the parameter e from zero, is it possible that the dis- 
cussed NNM ([7^ loses stability because of transforma- 
tion into a two-dimensional bush 



v{t)cN/4 + /i(t)c3^/4 



(76) 
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TABLE III: Asymptotic behavior of tlie critical parameters of 
NNMs in the thermodynamic hmit for the FPU-/3 chain with 
/3 > 0. 



NNM 



B[a^ j] 
_B[a*, ai] 

B[a^, iu] 



A, 



V 'iP 

■/2-K 

2tt 



ErJN 



3/3 JV^ 

3/3 JV^ 

27r 
3/3 JV 

3/3 JVi! 

27r 
3(3 JV^ 

27r 
3/3 JV2 



Refs. 



[I>i4-21J 

[3,1 
[3,2^ 

[2] 



TABLE IV: Asymptotic behavior of tire critical parameters of 
NNMs in the thermodynamic hmit for the FPU-/? chain with 
/3 < 0. 



NNM 


A, 




Refs. 


B[a^^ ai] 

B[a'^, m] 
B[a^ ^ai, a^u] 


1 2tt 

V 

V2n 


2n 


[2] 


3|/9|JV 

Z77 





3|/9|JV2 





A. The case /3 > 

Firstly, let us pay attention to the following interesting 
fact: four of six NNMs in the FPU-/3 chain, namely, [U 



B\a^,i] 



B[a3, tu] = v{t) ■ ^{1, ~2, 1 1 1, -2, 1 1 . . .}, 

B[a^u,ai] = • ^^7={o,i,i,o,-i,-i | 0,1,1,0,-1,-1 1 ...}, 

possess identical scaling of the stability threshold in the 
limit N 00: 



27r2 



(78) 



The scaling of Sc{N) for the 7r-mode 

B[a2,z]=Ki)--^{l,-l|l,-l| ■•■} (79) 



is exactly twice less than that determined by Eq. (|78p . 
Only for the mode 



B[a'^,ai] = iy{t) 



= {0,1,0,-1|0,1,0,-1| ...} (80) 



with two different functions of time vit), /i(t)? In gen- 
eral, such a bush describes not periodic, but a quasiperi- 
odic dynamical regime in the FPU-/3 chain. In contrast 
to the previous case, this stability loss scenario doesn't 
imply any extension of the vibrational modes set, but it 
means breaking the correlation ([75)) between two LNMs, 

i^N/4{t)^N/4 ^nd i^3N/4{t)'^3N/4- We don't discuss this 
scenario of the stability loss in the present paper. How- 
ever, our analysis allows us to assert that such scenario 
is inessential for all six symmetry-determined NNMs in 
the FPU-/3 chain. 



V. RESULTS AND DISCUSSION 

With the aid of the above discussed method, we inves- 
tigated scaling of stability thresholds in the thermody- 
namic limit — > 00 for all possible in the FPU-/3 chain 
nonlinear normal modes for both cases /? > and f3 < 0. 
The summary results are presented in Tables Hill II VI Let 
us comment on these results. 

In the first column of Tables IIIII IIVI the symbols of 
NNMs are given. In two next columns, we present the 
scaling laws for — 00 of stability thresholds for the 
amplitude of each NNM, Ac{N), and for its specific en- 
ergy, ec{N) (the energy per one particle of the chain). In 
the last column of the above tables, we give references 
to the papers in which stability thresholds of the corre- 
sponding NNMs are also discussed. 



the scaling law of the stability threshold turns out to be 
cardinally different: 



2lT 

3/3N' 



(81) 



A qualitative difference in scaling between this mode 
and all other NNMs for the FPU-/3 chain with /3 > can 
be seen in the left column of Fig. [TJ 



B. The case /3 < 

The stability properties of the same NNMs in the FPU- 
/3 chain with /3 < prove to be essentially different, as 
one can see in the right column of Fig. [T] 

Firstly, for three NNMs, namely, B[a2,z], B[a^,i] and 
B[a^, iu] the stability thresholds Sc{N) don't tend to zero 
when N 00. These modes lose their stability for a 
certain finite value Ac of the NNM's amplitude. For the 
above listed nonlinear normal modes, these values are 
equal respectively to 0.213, 0.092, and 0.138. 

The fundamental difference between scaling of edN) 
for the modes B[a'*,ai] and B[a'*,m] takes place in the 
case /3 < 0, as well as for the above discussed case /3 > 0: 
ec{N) - 1/iV for B[a^,ai] and - l/N^ for B[a'^,iu] 
(see Table ilVT). 

Studying of the stability threshold for the NNM 
B[a'^ii, ai] proves to be more difficult. In this case, the loss 
of stability is determined by the second scenario discussed 
at the beginning of the present section. Normalizing the 
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variational equations, described dynamics of the vibra- 
tional modes with basis vectors ^n/2i "^(^Af/6i^5Af/6)i 
simultaneously with the Duffing equation, corresponding 
to NNM B[a'^u,ai], we infer that for /3 < there exists 
an exponential detuning between the above modes for 
an arbitrary small amplitude of the investigated NNM. 
This means that ec{N) turns out to be equal to zero as 
indicated in Table HVl 

VI. SUMMARY 

In the present paper, a certain asymptotic technique 
for studying the stability loss of nonlinear normal modes 
in the FPU-/3 chain in the thermodynamic limit N ^ oo 
is developed. Using this technique we were able to obtain 



the scaling laws of the stability threshold Sc{N) for all 
six symmetry-determined NNMs that are possible in the 
FPU-/3 chain for both cases (3 > and /3 < 0. 

The general method [2^ for splitting the variational 
system for a given dynamical regime in a physical system 
with discrete symmetry into independent subsystems of 
small dimensions was applied for investigation of stability 
of NNMs in the FPU-/3 chain. The above dimensions for 
the considered case turn out to be equal to 1, 2 and 
3. This splitting allows us to construct numerically the 
stability diagrams (Fig. [!} that can help to reveal many 
interesting properties of nonlinear normal modes, such 
as qualitative behaviour of stability thresholds Sc{N) in 
the thermodynamic limit N ^ oo, the existence of the 
"second stability threshold" for some NNMs, existence of 
finite limits of edN) for certain modes, etc. 
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